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Abstract 

The Liouville equation, first Bogoliubov hierarchy and Vlasov equations with deriva- 
tives of non-integer order are derived. Liouville equation with fractional derivatives is 
obtained from the conservation of probability in a fractional volume element. This 
equation is used to obtain Bogoliubov hierarchy and fractional kinetic equations with 
fractional derivatives. Statistical mechanics of fractional generalization of the Hamil- 
tonian systems is discussed. Fractional kinetic equation for the system of charged 
particles are considered. 

PACS: 05.20.-y; 05.20.Dd; 45.10.Hj 

1 Introduction 

Fractional equations [1] are equations that contain derivatives of non- integer order [2l [3]. 
The theory of derivatives of non-integer order goes back to Leibniz, Liouville, Riemann, and 
Letnikov [3]. Derivatives and integrals of fractional order have found many applications in 
recent studies in mechanics and physics. In a short period of time the list of applications have 
become long. For example, it includes chaotic dynamics [H |5] , mechanics of fractal media 
P m [8], physical kinetics [H El [lOl [IH [12], plasma physics [H [H [Ej, astrophysics [16], 
long-range dissipation [T71 [TSl [12] , mechanics of non-Hamiltonian systems [201 [21] , theory of 
long-range interaction [221 [231 [21], and many others physical topics. 

In this paper, we derive Liouville equation with fractional derivatives with respect to co- 
ordinates and momenta. To derive fractional Liouville equation, we consider the conservation 
of probability to find a system in the fractional differential volume element. Using the frac- 
tional Liouville equation, we derive the fractional generalization of the Bogoliubov hierarchy 
equations. These equations can be used to derive fractional kinetic equations [H [H [TDl [12] • 
A linear fractional kinetic equation for the system of charged particles is suggested. 



In Sec. 2, we derive the Liouville equation with fractional derivatives from the conser- 
vation of probabihty to find a system in the fractional volume element of the phase space. 
In Sec. 3, we obtain the first Bogoliubov hierarchy equation with fractional derivatives in 
the phase space from the fractional Liouville equation. In Sec. 4, the Vlasov equation with 
fractional derivatives in phase space is considered. In Sec. 5, a linear fractional kinetic 
equation for the system of charged particles is suggested. Finally, a short conclusion is given 
in Sec. 6. 



2 Liouville equation with fractional derivatives 

A basic principle of statistical mechanics is the conservation of probability. The Liouville 
equation is an expression of this basic principle in a convenient form for the analysis. In this 
section, we derive the Liouville equation with fractional derivatives from the conservation of 
probability in a fractional volume element. 

For the phase space M^" with coordinates {x^, = (gi, ...,Pn), we consider 

a fractional differential volume element 

d'^V = d''Xi...d''X2n- (1) 

Here, is a fractional differential For the function f{x), 

2n 

dyix) = j2D:j{x){dxkr, (2) 
fe=i 

where D^^ is a fractional derivative [3] of order a with respect to Xk- The Caputo derivative 
[6] is defined by 

1 r 



D^f{x) = — / , ' ,\ dz, (3) 

where m — 1 < a < m, f^^\T) = d"^f{T)/dT^, and T{z) is the Euler gamma-function. For 
Caputo and Riesz Eq. [3J fractional derivatives, we have = 0, and D^^xi = {k ^ I). 
Using ([2]), we obtain 

d'^Xk = D^^XkidxkT, (a > 0). (4) 

Then 

{dx^r = {Dlxk)'' d-Xk. (5) 

For Caputo derivatives. 
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we get 

{dxkT = {D:,Xk)~' d'^Xk = r(2 - aX-^^Xfc. (7) 

The conservation of probability to find a many-particle system in the phase-space volume 
element d°'V may be expressed as 

- ^^^^^ = ^[/'(^' ^) ^S)] (8) 
for the usual volume element (a = 1), and 

_rf"V^^^^ = (i°[p(t,:r)(u,rf"S)] (9) 

for non-integer a. Here, p = p{t, x) is the density of probability to find a many-particle 
system in the phase-space volume element; u = u{t,x) is the velocity vector field, (i"S is a 
surface element, and the brackets ( , ) is a scalar product of vectors: 

2n 2n 2n 

u = J2^kek, ci"S = 5^d"5fce,, (u, ci"S) = ^ ^^^"5^, (10) 

k=l k=l k=l 

where are the basic vectors of Cartesian coordinate system, and 

d'^Sk = d°'xi...d°'Xk-id'^Xk+i---d°'x2n, (11) 

The functions Uk = Uk{t, x) define Xk components of the velocity vector field u{t, x), the rate 
at which probability density is transported through the area element ^"5*^. In the usual case 
[a = 1), the outflow of probability in the Xk direction is defined by 

d[puk]dSk = D^^[puk]dxkdSk = D^^[pUk]dV. (12) 

The fractional generalization of Eq. (fT2|) is 

c?"[pMfc]rf"5fc = D:^[puk]{dxrd^Sk = 

= (Dlxk)-' Dl[puk]d^Xkd'^Sk = {Dlxk)-' Dl[pUk]d'^V. (13) 
Here, we use Eqs. f lTT]) . ([T]) and (jS]). Substitution of Eq. f[T^ into Eq. gives 

o 2?! 

- d'^V^ = d-Vj2{Dlxk)-'Dl[puk]. (14) 



dt 

As a result, we obtain 

dp 



k=l 



2n 



k=l 



For Caputo derivatives 

dt 



o 2n 

p -n2-a)Y,^t'mpM- (16) 



k=l 

Equation (fTSi) is the Liouville equation that contains the derivatives of fractional order a. 
Fractional Liouville equation ( |T5l) describes the probability conservation to find a system in 
the fractional volume element ([T]) of the phase space. 

For the coordinates = (gi, g„,pi, Eq. (fT5|) is 



rj n n 



fc=l ' ■ k=l 



where Vk = Uk, and Fk = Uk+n {k = l,...,n). The functions Vk = Vk(t,q,p) are the 
components of velocity field, and Fk = Fk{t,q,p) are the components of force field. In 
general, 

D^^ [pFk] ^ pD^Fk + FkD<^^^p. (18) 

Suppose that Fk does not depend on pk, and the kth. component Vk of the velocity field does 
not depend on fcth component qk of coordinates. In this case, Eq. f|T71) gives 



dt 



r\ II It 

" + 11^^ {Dl<lk) DIP + J2Fk {D-^Pu) D^p = 0. (19) 

fc=l k=l 



If we consider the fractional generalization of Hamiltonian system [21] , then Vk and Fk can 
be represented as fractional derivatives of some function H{q,p): 

Vk = D^^Hiq,p), Fk = -DlH{q,p). (20) 

For a = 1, we have the usual Hamiltonian system. 

The classical system that is defined by the equations 

^ = Vk{t,q,p), ^ = Fk{t,q,p) (21) 

is called a Hamiltonian system, if the right-hand sides of equations (12T|) satisfy the Helmholtz 
conditions 

dV dV dV dFi dFi dF. , , 

i I = 0, ^ H - = 0, - - = 0. (22) 

dpj dpi ' dqi dpj ' dqj dqi 

In this case, Eq. f l2T]) can be presented in the form 

dqk dH dpk dH 
dt dpk dt dqk 



(23) 



that is uniquely defined by the Hamiltonian H. The fractional Hamiltonian systems are 
defined in Ref. pO]. 

For Eq. ([20]), equation is 

n n 

Yt+Y. iDl^k) D;HDIp - J2 {D;,Pk) DlHD-^p = 0. (24) 

k=l k=l 

We can define 

n 

{A,B}^ = J2 {{Dlq,)-' DIAD'^^B - {D^Pky' D^BD^a) . (25) 

k=l 

For a = 1, Eq. (125|1 defines Poisson brackets. Note that the brackets (l25l) satisfy the relations 

{AB}^ = -{B,A}^, {1,A}„ = 0. 

In general, the Jacoby identity cannot be satisfied. The property {l,y4}Q, = is satisfied 
only for Caputo and Riesz fractional derivatives (-D"l = 0). For the Riemann-Liouville 
derivative, D^l ^ 0. Using Eq. fl2^ . we get Eq. fl2^ in the form 

1^ + {p, = 0. (26) 

As a result, we have the Liouville equation for the fractional generalization of Hamiltonian 
systems |2T] that are defined by Eq. (!20|) . For a = 1, Eq. (126|) is the usual Liouville equation 
for Hamiltonian systems. 



3 Bogoliubov equation with fractional derivatives 

Let us consider a classical system with a fixed number of identical particles. Suppose 
the /cth particle is described by the generalized coordinates g^s and generalized momenta 
Pfcs, where s = l,...,m. We use the notations = {qki, ■■■,qkm) and = {pki, ■■■,Pkm)- 
The state of this system is described by the N-particle density of probability pn in the 
2mA^-dimensional phase space: 

P7v(q, P, t) = p(qi, pi, Qat, Pat, t). 

The fractional Liouville equation is 



dpN _ ^ 
dt 

k=l 



5^(d^JV,P^) + D°jF,p^)), (27) 
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where 

m 



m 



Dp.F. = {D-,,P,r' D-^F, = ^ (Dl^P^s)-' D^F,,. (29) 

s=l 

The one-particle reduced density of probabihty pi can be defined by 

Pi(q,P,^) =p(qi,Pi,^) = I[2,...,N]pN{q,p,t), (30) 

where J[2, ...,N] is an integration with respect to variables q2, qAr, P2, Pn- Obviously, 
that one-particle density of probability satisfies the normalization condition 

/[l]pi(q,p,t) = l. (31) 

The Bogoliubov hierarchy equations [271 1211 [291 EO] describe the evolution of the reduced 
density of probability. They can be derived from the Liouville equation. Let us derive the 
first Bogoliubov equation with fractional derivatives from the fractional Liouville equation 
(127]) . Differentiation of Eq. (I30l) with respect to time gives 

Using Eq. fl271) . we get 

N 



dpi 
dt 



-/[2, ...,iV] ^(d^JV.p^v) + D^jF.p^v)). (32) 



k=l 



Let us consider the integration over and for k-particle term of Eq. (!32l) . Since the 
coordinates and momenta are independent variables, we derive 

^ _1 / \ +00 

/[q,]D°jV,p;v) = I[qk] {D^q,)-' D"^^{y,pr,) ~ V,p^ = 0. (33) 

\ / — oo 

For example, the Caputo derivatives give 

/[qfc]D^jVfep^) = r(2 - a)i[q,]cit-'D^^iV,p^) = 

= r(a)r(2 - a)/"[qfc]D-jVfcp^) = r(a)r(2 - a) (v.pjv)^" = 0, (34) 

V / — oo 

where /"[q^] is a fractional integration with respect to variables q^. In Eq. (1331) . we use the 
fact that the density of probability p^v in the limit q^ — > ±oo is equal to zero. It follows 
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from the normalization condition. If the hmit is not equal to zero, then the integration over 
the phase space is equal to infinity. Similarly, we obtain 



Then all terms in Eq. fl521) with k = 2, ...,N are equal to zero. We have only the term with 
k = 1. Therefore Eq. (132!) has the form 

^ = -/[2,...,iV](D^^(Vip^) + D^^(Fip^)). (35) 

Since the variable qi is an independent of q2, Qat and p2, Pat, the first term in Eq. 
fl35l) can be written as 



J[2, iV]D^jVipAr) = D^^ Vi/[2, iV]p^ = D^^(Vipi)- 
The force Fi acts on the first particle. It is a sum of the internal forces 

Fife = F(qi,pi,qfc,pfc,t), 
and the external force F^ = F^(qi, pi, t). In the case of binary interaction, we have 

N 



Fi = F^ + 5^Fifc. (36) 



k=2 



Using Eq. (1361) . the second term in Eq. (!35|) can be rewritten in the form 



N 



/[2, iV]D^^ (Fip^) = i[2,...,N] (d^^ (F^p^) + ^p. (FifcP^) 

k=2 

N 

= D^^ (F^pi) + ^ B-J[2, N] (Fifep^) . (37) 

k=2 

We assume that pn is invariant under the permutations of identical particles [3T] : 

PAr(...,qfc,Pfc, ...,q«,p«, ...,t) = pN{---,qi,Pi, ...,qk,Pk, ■■■,t). 

In this case, pa? is a symmetric function, and all (A^ — 1) terms of sum fl371) are identical. 
Therefore the sum can be replaced by one term with the multiplier [N — 1): 

N 

^/[2,...,iV] D^^^(Fifcp^) =(iV-l)/[2,...,iV] D^^(Fi2Piv). (38) 



Using /[2, A^] = /[2]/[3, A^], we rewrite the right-hand side of (l38l) in the form 

i[2] D^^ (Fi2/[3, iV]p^) =D^^/[2] (F12P2), (39) 

where p2 is two-particle density of probabihty that is defined by the fractional integration of 
the iV-particle density of probability over all and p^, except k = 1,2: 

P2 = p(qi, Pi, q2, P2, t) = i[3, N]pN{q, p, t). (40) 

Since pi is independent of q2, P2, we can change the order of the integrations and the 
differentiations: 

J[2] D^^(Fi2P2) =D°^/[2]Fi2P2. 



Finally, we obtain 

dp 



i-t-D^^(Vipi) + D^^(F^pO =/(p2). (41) 



dt 

Here, /(P2) is defined by 

/(p2) = -(iV-l)D^^/[2]Fi2P2. (42) 

The term /(P2) describes a velocity of particle number change in 4m-dimensional two-particle 
elementary phase volume. This change is caused by the interactions between particles. 

Equation (14T|) is the fractional generalization of the first Bogoliubov equation. If a = 1, 
then we have the first Bogoliubov equation for non-Hamiltonian systems [36] . For Hamilto- 
nian systems, 

p ^ dH{qi,pi) ^ dH{qi,pi) 

dqi ' dpi 
and Eq. dH]) has the well-known form [271 EHl [29l ED] . 

4 Vlasov equation with fractional derivatives 

Let us consider the particles as statistical independent systems [32l[33]. Then 

P2(qi,Pi,q2,P2,t) = pi(qi,pi,t)pi(q2,P2,t)- (44) 
Substitution of Eq. (jS]) into Eq. gives 

/(P2) = -D-^piJ[2]Fi2Pi(q2,P2,t), (45) 
where pi = pi(qi, pi, t). As a result, the effective force is 

F^^^(qi,pi,t) =/[2]Fi2Pi(q2,P2,t). 
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In this case, we can rewrite Eq. ( H5l) in the form 



/(P2) = -D^,(piF^^0- (46) 
Substituting of Eq. fH6l) into Eq. (j4Tl) . we obtain 

^ + D^,(Vipi) + D^^ ((F^ + (AT - 1)F^^0pi) = 0. (47) 

This equation is a closed equation for one-particle density of probability with the external 
force and the effective force F^'^-^. Equation ( l47l) is the fractional generalization of the 
Vlasov equation that has coordinate derivatives of non-integer order. For a = 1, we get the 
Vlasov equation for the non-Hamiltonian systems that is described by non-potential fields. 
For Hamiltonian systems (H3|) . equation (HTl) has the usual form [32| [33]. 



5 Linear fractional kinetic equation for a system of 
charged particles 

Let us consider fractional kinetic equation (HT]) with /(P2) = 0, Vi = p/m = v, and F"^ = eE, 
B = 0. In this case, Eq. fHTl) has the form 



df 
dt 



^ + (v,D^/)+e(E,D^/) = 0, (4J 



where / = pi is one-particle density of probability, and 

m 

(v,D^/) = 5^(^„D°J). (49) 



If we take into account the magnetic field (B 7^ 0), then we must use the fractional 
generalization of Leibnitz rules: 

Dg(/«) = E r(. + i)r(^-. + i) (°r7)i^;«. (50) 

where r are integer numbers. In this case, Eq. fHS]) has the addition term 

— D^([p,B]/) = ^ VD^^ {SMrnPlB^f) = — y2eMrnBmr>;^ (pif) = 

mc ^ mc ^-^ ' mc ^-^ ^ 

klm klm 
= y^l^klmBm — [ ^ ^ — — [Dp"V]4iP[ = 

mc ^-^ i r + 1 i a — r + 1 

klm r=0 \ / \ I 



e 

mc 



klm 

'_Y.^Hr.Br^Pl[^;j] = ^ mj), [P, B]) . (51) 

klm 

Let us consider the perturbation [SU [35] of the density of probabihty /o in the form 

f{t,q,p) = fo + 5f{t,q,p), (52) 

where /o is a homogeneous stationary density of probabihty that satisfies Eq. PHI) for E = 0. 
Substituting of ( |52l) into Eq. ( HHl) . we get for the first perturbation 



^ + (v, D^5/) + e(E, D^/o) = 0. (53) 

Equation fl53l) is the hnear fractional kinetic equation for the first perturbation of the density 
of probabihty. Solutions of fractional linear kinetic equations of type (153|) are considered in 
Ref. [10]. For E = 0, the function 6f is described by 



(^7.t)-i/°L„ [qsigst)-'/''] , (54) 

where gs = Vs{D^^qs)^^ , and 

1 r+oo 

L^[x] = — / dfc e-'^'e-l*^!'' (55) 



is the Levy stable density of probability [40] . 

For a = 1, the function fl55|) gives the Cauchy distribution 



TT + 1 ' 

and Eq. is 

1 (^7.^)-' 



For a = 2, we get the Gauss distribution: 

1 



(57) 



and the function (IMl) is 



L,[x] = ^e-^^/^ (58) 

2a/ TT 



(^3t)-V2j_e-^^^/(^^»*). (59) 
2Vvr 
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For 1 < a < 2, the function La[x] can be presented as the expansion 



LJx] = > (-x) — H — ^sin(n7r/2). 

n=l 



(60) 



The asymptotic (a; — oo, 1 < a < 2) is given by 



1 



oo 



r(l + na) 
n! 




n 



X 



—na 



sin(n7r/2). 



(61) 



n=l 



As a result, we arrive at the asymptotic of the solution, which exhibits power-like tails for 
X — oo. The tails is the important property of the solutions of fractional equations. 

6 Conclusion 

In this paper, we consider equations with derivatives of non-integer order that can be used 
in statistical mechanics and kinetic theory. We derive the Liouville, Bogoliubov and Vlasov 
equations with fractional derivatives with respect to coordinates and momenta. To derive the 
fractional Liouville equation, we consider the conservation of probability to find a system in 
the fractional differential volume element. Using the fractional Liouville equation, we obtain 
the fractional generalization of the Bogoliubov hierarchy equations. Fractional Bogoliubov 
equations can be used to derive fractional kinetic equations [H [9l [10]. The fractional kinetic 
is related to equations that have derivatives of non-integer orders. Fractional equations 
appear in the description of chaotic dynamics and fractal media. For the fractional linear 
oscillator, the physical meaning of the derivative of order a < 2 is dissipation. Note that 
fractional derivatives with respect to coordinates can be connected with long-range power- 
law interaction of the systems [22l [231 [21] ■ 
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